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ABSTRACT 

Super-Eddington accretion onto black holes (BHs) may occur at Ultra-Luminous compact X-ray sources in 
nearby galaxies, Galactic microquasars and narrow-line Seyfert 1 galaxies (NLSls). Effects of electron scattering 
(opacity and Comptonization) and the relativistic correction (gravitational redshift and transverse Doppler effect) 
on the emergent spectra from super-Eddington accretion flows onto non-rotating BHs are examined for 10'^ and 
masses (Mbh)- With in [= M/(LEdd/c^), where M is the accretion rate] > 100, the spectral hardening 
factor via electron scattering is ^ 2.3-6.5. Due to the m-sensitive hardening factor, the color temperature of the 
innermost radiation is not proportional to L"^^, differing from the simplest standard accretion disk. The model 
is applied to optical-soft X-ray emission from NLSls. We pick up one NLSl, namely PG 1448-1-273 with an 
inferred Mbh of 10'''*Mq, among the highest m candidates. The broadband spectral distribution is successfully 
reproduced by the model with an extremely high m (=1000) and the viscosity parameter a of 0.01. This implies 
that this object, as well as some other highest rh systems, is really young: the inferred age, Mbh/M, is about 10^ 
years. We also briefly discuss the distribution of m for transient and highly variable NLSls, finding that those are 
located at 3 ^ m ^ 300. Such a moderately high accretion rate is indicative of thermal instability. Furthermore, m 
for a possible type-2 counterpart of NLSls, NGC 1068, is found to be similar to m for NLSls. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — galaxies: nuclei — 
X-rays: galaxies — X-rays: stars 



1. INTRODUCTION 

Super-Eddington accretion phenomena onto black holes 
(BHs) may be important for formation and evolution of su- 
permassive BHs. It may have been much more common in 
the young universe, when the average black hole mass (Mbh) 
would have been smaller than now. Such phenomena are possi- 
ble sources of UV-Extreme UV (EUV) photons that may have 
governed the thermal and star formation history of the universe. 
Moreover, those extreme accretion rate phases provide us a 
laboratory in examining our knowledge on accretion phenom- 
ena, which have been studied intensively under sub-Eddington 
regimes (e.g., Shakura & Sunyaev 1973). 

Some Ultra-Luminous compact X-ray sources (ULXs) and 
Galactic microquasars show large bolometric luminosity ( ^ 
the Eddington luminosity LEdd of 1.4 Mq) and high color tem- 
perature (~ 1 keV), often referred as Tjn (e.g., Roberts & War- 
wick 2000; Makishima et al. 2000; Colbert & Ptak 2002). We 
note that recent results of ULXs obtained with XMM-Newton 
seem to show lower Tjn (^--^0.15 keV; Miller et al. 2003). Some 
possible physical explanations for their nature, including rotat- 
ing BHs, high accretion rate (M > Lndd/c^, where M is the gas 
accretion rate), intermediate BH mass (^lOOM©) and beamed 
radiation (see, however, Misra & Sriram 2003) have been dis- 
cussed. 

There is a subclass of Active Galactic Nuclei (AGNs) that 
are candidates for super-Eddington accretion, namely narrow- 
line Seyfert 1 galaxies (NLSls) and their luminous counter- 
part, narrow-line QSOs. They are characterized as follows (see 
Pogge 2000). (i) They have narrower Balmer lines of hydrogen 
(e.g., FWHM of H/? < 2000 km s"'), relative to usual broad- 



line Seyfert 1 galaxies (BLSls) and QSOs having FWHM of 
H/3 > 5000kms-i (e.g., Osterbrock& Pogge 1985). (ii) They 
often emit strong optical Fe II multiplets (e.g., Halpern & Oke 
1987). (iii) Their optical-soft X-ray, big blue bump is hotter 
than BLSls/QSOs: they have bluer optical spectra (Grupe et 
al. 1998) and often show steep and luminous soft X-ray ex- 
cess (Pounds et al. 1996; Otani et al. 1996; Boiler et al. 
1996; Wang et al. 1996; Laor et al. 1997; Leighly 1999b). 
(iv) Rapid soft/hard X-ray variability is another characteristic 
of NLSls (Otani etal. 1996; Boiler et al. 1997; Leighly 1999a; 
Hayashida 2000). The soft X-ray spectral variability of NLSls 
is opposite to that of BLSls: spectra get harder when they be- 
come brighter (Cheng, Wei & Zhao 2002). 

When the luminosity from the accretion disk/flow L becomes 
close to Lfidd L3 x 10'^^(MBH/10^MQ)erg s"'], it has been 
shown that the advective energy transport in the flow dominates 
over radiative cooling (Abramowicz et al. 1988; section 10.3 of 
Kato, Fukue & Mineshige 1998). Such a disk is called an the 
optically-thick ADAF (advection-dominated accretion flow) or 
slim disk, since it is moderately thick geometrically. Abramow- 
icz et al. showed that L saturates at around a few times Ledd 
when M exceeds Lndd/c^, and that the flow shines even in- 
side the marginally stable orbit [3rsch for non-rotating BHs; 
where rsch is the Schwarzschild radius and rsch = 2GMbh/c^ ~ 
3 X 1O'2(Mbh/1O^M0) cm = 10"^ (Mbh/IO^Mq) It-days.]. The 
physical reason why the advection turns on at super-Eddington 
regime will be the fact that photons can not escape from the 
flow before photons/gas are swallowed by the central BH (the 
photon trapping effect; Begelman & Meier 1982; we discuss it 
in §6.2). Slim disk models for ULXs have been constructed to 
explain High Tin (Watarai et al. 2000) and mysterious behaviors 
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in Tin v.s. Lx diagrams (Watarai, Mizuno & Minesliige 2001). 
Similarly, these models are successful in discussing strong soft 
X-ray emission (large Tin relative to other AGNs) and time vari- 
ability of NLSls (Mineshige et al. 2000, with Mbh and M be- 
ing free parameters). Currently, the most promising picture of 
NLS 1 s and narrow-hne QSOs is that they contain relatively less 
massive black holes (with Mbh ~ lO^^^M©) and higher L/L^a 
(with L ~ LEdd)- 

As to the emergent spectra from super-Eddington accretion, 
Szuszkiewicz, Malkan & Abramowicz (1996) and Wang et al. 
(1999) considered the opacity due to electron scattering, but 
without Comptonization nor the relativistic correction. Mi- 
neshige et al. (2000) and Watarai et al. (2000) simply apply the 
local blackbody approximation, whereas Watarai et al. (2001) 
take the spectral hardening factor of 1 .7 from Shimura & Taka- 
hara (1995) in order to mimic the effect of Comptonization. We 
note that this factor is derived in the sub-Eddington regime, and 
thus the extent of Comptonization at super-Eddington phase is 
not evaluated by it (see §3.2). Self-irradiation and -occultation 
are discussed by Fukue (2000) with self-similar solutions of 
the slim disk model (e.g., Begebnan & Meier 1982). Most re- 
cently, Wang & Netzer (2003) have studied the effect of Comp- 
tonization with self-similar solutions, in order to explain X-ray 
emission from NLSls. 

A big problem remains when we compare the slim disk 
model with the observed optical/UV-soft X-ray spectra of 
NLSls. Namely the predicted soft X-ray is too step compared 
with observations (Szuszkiewicz et al. 1996; Wang et al. 1999). 

We examine the effects of electron scattering in opacity, of 
Comptonization and of the relativistic correction upon spectra 
from transonic slim accretion flows onto non-rotating BHs, in 
order to figure out to what extent those effects change the emer- 
gent spectrum and whether the effects are enough to fit the ob- 
served spectral energy distributions. Now that Mbh estimations 
for AGNs are accumulating, we are able to compare the models 
and observational data with less free parameters than previous 
studies. We have obtained a promising fit to the broad-band en- 
ergy distribution for the first time. Methods and assumptions 
of the numerical calculations are given in §2. We then discuss 
the physical quantities and emergent spectra of the slim disk 
in §3. We turn to the comparison of the model with observed, 
broad-band spectra of NLSls in the next section. The growth 
timescale of supermassive BHs is discussed in §5. In §6, sev- 
eral discussions are presented. The final section is devoted to a 
summary. Throughout the present study we use the normaliza- 
tion MEdd(= ^Edd/ c^) for M: 

MEdd ^ 1.3 X 10'^Mbh/IOMq) g S-' 
~ 1.3 X 10^4 (Mbh/ lO'M0)gs-' 
~ 2.0 X 10-2(Mbh/10''Mq) Mq yr-\ (1) 

Hereafter, m refers to M/MEdd- 

Throughout this study, the luminosity is calculated from ob- 
served flux assuming isotropic radiation, zero cosmological 
constant, deceleration parameter = 0.5, and Hubble constant 
Hq = 75 km/s/Mpc. 

2. NUMERICAL PROCEDURES 

2. 1 . The Structure of the Accretion Disk 

In order to discuss the shape of the emergent spectra from 
the iimermost regions, it is essential to solve the full, derivative 
equations, rather than the self-similar solutions or other approx- 
imate solutions (based on a priori iimer boundary conditions or 



no regularity conditions). The transonic nature of the accretion 
flow should also be carefully treated. 

We solve the steady-state, transonic accretion disk structure. 
The numerical methods are basically the same as those adopted 
in Matsumoto et al. (1984). The code has been extended for 
time-dependent simulations (Matsumoto, Kato & Honma 1989; 
Honma et al. 1991a,b) and it has been made use of in subse- 
quent papers; Takeuchi (2000), Watarai et al. (2000, 2001) and 
Mineshige et al. (2000). Here, we summarize the numerical 
procedures, and briefly describe the basic equations in the fol- 
lowing paragraphs. We employ the pseudo-Newtonian poten- 
tial, V =-GM/(r- rsch) (Paczynski & Witta 1980), and cyUn- 
drical coordinates, (r,(p,z)- 

In this study, we use vertically integrated (height-averaged) 
derivative equations. The scale-height of the disk H is compa- 
rable to radius r in the case of extremely large accretion rates 
M ^ lOOLEdd/c^. One may wonder if such a height-averaged 
approach would not be appropriate. The same problem arises 
in an optically thin ADAF (Narayan & Yi 1995; see section 
10.2 of Kato et al. 1998 and references therein). Narayan 
& Yi (1995) found that the 2-dimensional solutions of the ex- 
act non-height-integrated equations agree quite well with those 
of the simplified height-integrated equations when the 'height- 
integration' is done along a constant spherical radius, rather 
than along z at constant cylindrical radius (see Narayan 1997). 
The height-integrated equations therefore are fairly accurate de- 
scriptions of quasi-spherical, advective flows. 

We construct the vertically integrated equations, using the 
integrated variables, such as surface density, S = J pdz = 2pH, 
and integrated total (gas plus radiation) pressure, H = J pdz = 
2pH. Throughout this study (except for the absorption coef- 
ficient k), each quantity with a bar on top means the vertically 
averaged value, evaluated from the value at the mid plane (value 
with the suffix "mid"): p = (16/35) pmid, P = (128/315)/7mid 
and r = I Tjnid, respectively (Hoshi 1977; see Matsumoto et al. 
1984). The scale-height of the flow H is determined by hy- 
drostatic balance between the vertical gravity and pressure: i.e. 
H = Ok/c„ where is the sonic speed. Equations for the con- 
servations of the mass, momentum, and angular-momentum are 

-IwrJ^Vr = M = const. , (2) 

dvr idn_ f-el n dlnQK 

and 

m-lin) = -2TTr^Tr^, (4) 

respectively. Here, Q.{=v^/r), 9.Yi{= {GM /rfl^ /{r-rs^^)], £ 
(= fv^), (= '■^^^k) and £ia are respectively, the angular fre- 
quency of the gas flow, the Keplerian angular frequency in the 
pseudo-Newtonian potential, the specific angular momentum, 
the Keplerian angular momentum and the specific angular mo- 
mentum finally swallowed by the BH. As to the viscous stress 
tensor we adopt the usual prescription Tr^ = -all with the vis- 
cosity parameter a. Unless Otherwise noted, we apply a of 0.1. 
We neglect self-gravity of the disk, for simplicity. We will see 
that this is self-consistent with the computed density profile. 

The energy budget (per unit surface area) at each radius is 
symboUcaUy written as 

Gadv = Gvis ~ 2rad> (5) 

where Q^^yloc llVrT{ds/dr)] is the advective cooling with s be- 
ing specific entropy, and the other two terms on the right-hand 
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side represent viscous heating and radiative coohng, respec- 
tively. Assuming that the effective optical depth Teff is larger 
than unity, we apply the usual, diffusion approximation for ^ 
(cf. §3.1); 



(6) 



where the total opacity coefficient is «; = Kabs + K^s- Here, Mrad 
and tdiff mean the radiation energy density per unit volume 
{af^) and the diffusion time scale of photons (see §3.1), re- 
spectively. This equation describes the photon leakage. Neu- 
trino leakage in proto-neutron stars is discussed in §6.2. If Teff 
gets smaller than unity, Mrad will be lower than aT'^. The cur- 
rent cooling rate (eq. 6) will then overestimate the true cooling 
rate (Beloborodov 1998). An improvement of the formula, e.g., 
multiplying the cooling rate by (1 -exp"^"^" )^ would be required 
in the future. 

The outer boundary conditions are imposed at r = 2.0 x 
lO^rsch, where each physical quantity is taken from the formula 
of the standard disk (Shakura & Sunyaev 1973). The choice of 
the outer radius does not influence on the inner solutions (e.g., 
Abramowicz et al. 1988). The basic equations are integrated by 
the semi-impUcit method from the outer boundary to the irmer 
one taken at r = 2.5rsch for m < 10 and 1.01 rsch for m > 100. 
The solutions satisfy the regularity condition at the transonic 
radius, at about 2.7rsch, and we adopt the free boundary condi- 
tions at the inner edge. 

The system is taken to be a disk consisting of thermal plasma 
around a Schwarzschild black hole of mass Mbh. We also as- 
sume 1 -Temperature gas (i.e., temperatures of protons and of 
electrons are the same). The Rosseland mean opacity is adopted 
for the absorption opacity Kabs; 'tabs = ^abs.o 

pT''^-^. In this 

study, every optical depth is measured between the surface and 
the mid-plane of the flow (not from one surface to the other). 
Generally, the optical depth of absorption Tabs (HpK^hs) is much 
less than that of electron scattering Tes (Hpn^s, with the elec- 
tron scattering opacity k^s of 0.4 cm^/g) in the region of inter- 
est (r< 10^"Vsch; cf. §3.2 in Kato et al. 1998). The flow is 
optically thick; Tabs + Tes 3> 1 ■ The absorption opacity is an im- 
portant factor when we are considering whether the effective 
optical depth of the flow [Teff = v^Sr^bsC^Vbs + r^; Rybicki & 
Lightman 1979] exceeds unity or not. To estimate the absorp- 
tion coefficient Kabs.o (mainly bound-free transitions for AGN 
disks with solar abundances; see e.g., Laor & Netzer 1989), 
we simply take 30 times the free-free absorption coefficient 
(«ff,o = 6.4 lO'^^ cm^/g) for the super massive BHs, as is done 
in Czemy & Elvis (1987). For much smaller Mbh (32 Mq), 
we adopt Kff.o as Kabs.o- The timescale of Comptonization and 
that of Coulomb colhsions are much shorter than the accretion 
timescale (Beloborodov 1998). 

We neglect the gas evaporation from the disk (e.g., Meyer & 
Meyer-Hofmeister 1994), since it becomes weak when the gas 
accretion rate approaches the Eddington limit (Liu et al. 2002). 
The coronae above and beyond the disk are not included in this 
study. It will be necessary to consider them when applying the 
models to hard X-ray spectra (Wang & Netzer 2003). 

In total, the input parameters required for the calculations are 
Mbh, M and a. 

2.2. Spectral Calculations 

We follow Czemy & Elvis (1987) and Wandel & Petrosian 
(1988), in dealing with the effect of opacity of electron scat- 
tering (i.e. the modified blackbody) and the effect of the 



energy exchange via Compton scattering (Comptonization). 
More precise treatments, including such as the effects of ra- 
diative transfer (see Shimura & Manmoto 2003), will be re- 
quired for more detailed comparison of models and observa- 
tions. Given the T^ff, the local spectrum /^(r) can be expressed 
by Uir) = Bv{T)fv{T,r), where B^iT) is the Planck function. 
Here, fv{T,r) is 

2|l-exp-^''''^'"*| 



MT,r) = 



(7) 



1 + \/[nhsAT, r) + Tes(r)] / T^hsAT, r) ' 

where Tabs,i/(r,r) is the optical thickness for absorption with a 
given temperature T at frequency v, and the effective optical 
depth Teff,,y(r,r) is computed as .J^T,,b^AT:r)[T^jB^,JT,r) + T^,(r)\. 
The parenthesis including the exponential term takes into ac- 
count the finite, effective optical depth of the flow. For low 
Teff,;y(r, r) (^ 1), it describes the thermal bremsstrahlung from 
an optically thin medium. If the absorption opacity dominates 
over the opacity of electron scattering (and if the flow is effec- 
tively, optically thick), f„(T,r) is unity. The color temperature 
T is determined so that the surface emissivity does not change: 
dvl^ir) = acT^fi/4n. 
A simple evaluation of the spectral hardening factor due to 
the modified blackbody is (Madej 1974): 



Tabs 



(8) 



, Tabs + Tes j 

Compton scattering is important when the Compton param- 
eter y = (Akf /niec^) Maxij-^^ .t^^) is greater than 1. However, below 
the last thermalization surface, where the effective optical depth 
measured from the surface equals unity, the energy exchange 
via Compton scattering will be statistically canceled due to the 
efficient absorption. Scattering above the surface is in question. 
Then, the relevant opacity for electron scattering, t^, will be 

Tes /Teff for Teff > 1 and Tes for Teff < 1: 

r' = — (9) 

Max(Teff,l) ^ ' 

Finally, the effective Compton-j parameter is 



^Max(Tes,Tes). 



(10) 



Comptonization will considerably modify the emergent spec- 
trum if > 1 . The radial distribution of increases inwards 
steeply, reaching a maximum at ^ (3 - 5)rsch, and it decreases 
inward (Wandel & Petrosian 1988). In the regime for which 

Tes > Tabs and Teff > 1, 

y^oir-^p —r^ , (11) 

V l^fCfi ) 

yielding a strong temperature dependency. 

Photons generated at large optical depth undergo many scat- 
terings and change their energy as a result of Comptonization. 
The effect on emergent spectra can be estimated by assuming 
that a fraction ffa ^ of the photons at frequency v are shifted and 
cumulated into a Wien peak of an average energy 3^r. For a 
given frequency v, the fraction of these thermalized photons to 
all generated photons /th,i/ is given as 

/th,!/ = 



exp 



(12) 



Te^ ln[ 1 + 4fcr/OT,c2 + 1 6(/fer/OT,c2)2] 

(Svensson 1984), if the absorption is negligible. For high- 
energy photons (fer < hv), /th,,^ would be able to exceed unity. 
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1 10 100 1000 lOM 10 100 1000 10" 
Radius (Rsch) Radius (R^ch) 

Fig. 1 . — Opacity of electron scattering (tcs, solid lines), that of absorption (r^bs. dotted lines) and the effective opacity (Teff, dot-dashed lines). Here, Mbh and 
a are chosen to be IO^^Mq and 0.1, respectively. It is shown that Tcs is larger than Tat, by several orders at r ;S lOrsch for all the m (1-1000). The left panel is 
drawn for ih =10 and 1000, while the right one for ni = 1 and 100. Open circles indicate the position where Prad = Pgas- Inner parts are p,ad -dominated regions. With 
m = 1000, entire region presented here are Prad"dominated. 



We therefore put a constraint so that fth,i, < 1 .0. Since the re- 
distribution of low-energy photons into Wien peak is the prob- 
lem, this artificial constraint on high-energy photons will be a 
reasonable approximation (cf. Wandel & Petrosian 1988). The 
formula (eq. 12) can be adapted to a more general case if Tes in 
the equation is replaced by the total optical depth r^"'"' (Czerny 
& Elvis 1987): 

total Tes + Tabs M 
1+Teff 

« Tes (for Teff < 1 , and Tes > Tabs) 

~ tU^t'II'^ (for Teff > 1 , and Tes > Tabs). (13) 

The removed photons will result in a total additional energy 
output (we use 2.1 kT instead of 3fer, following Wandel & Pet- 
rosian 1988) 

/,h,.=2.7^r/ dvf,^^,^. (14) 

The distribution of the shifted photons around the Wien peak 
is approximated by a blackbody spectrum with a normalization 
constant C given by 

lA,.yK.CacT^ (15) 

As a whole, the deviation of the emitted spectrum from a black- 
body is finally described as 

/.(r,r) = 



2 [ 1 - exp 



--TrffCr.r) 



1 + \/[T"abs,i'(7', r) + Tes(?-)] / Tabs,i/(7', r) 



[l-/th,.]+C, (16) 



instead of fy{T,r) in eq. (7). Again, the temperature T (> Teff) 
is adjusted to reproduce the original surface emissivity, so that 
the luminosity of the flow is conserved. As the result, the 
emergent spectra become harder than spectra based on the lo- 
cal blackbody approximation [B^(reff); Czerny & Elvis 1987; 
Wandel & Petrosian 1988]. A hot surface layer will arise whose 
temperature T is much larger than Teff, as seen in standard disks 
with m> \ (Shimura & Takahara 1993). 

Next, we add the effects of the relativistic correction on the 
effective temperature T^n and emergent spectra. Since only 



the gravitational redshift and transverse Doppler shift (i.e. no 
Doppler boosting) are considered, the computational results 
correspond to face-on views of the flow. We assume that the 
rotational velocity is almost the same as the Keplerian veloc- 
ity vk. The consistency will be shown later with the calculated 
velocity field (§ 3.1). We employ the relativistic correction for 
the Keplerian disk (e.g., Kato et al. 1998, §3.5.3); 

We take l+z„. = [Max(10-\ 1 - 3rsch/2^)]"'/^ for the mo- 
ment. The effective temperature Teff becomes lower by (1 + 
Zf.c )"' [i.e., emissivity per surface area decreases by (l+Zr.c.) "^], 
compared with the cases without the relativistic correction (e.g., 
Watarai et al. 2000, 2001; Mineshige et al. 2000). 

We have applied this expression to describe the shape of the 
radiation spectrum emitted at every radius of the disk. 

3. NUMERICAL RESULTS 

In this section, we show the basic physical quantities in the 
accretion flow and emergent spectra for Mbh of lO^^M© and 
a = 0.1. For AGNs (Mbh ^ IO^Mq), Mbh is now estimated by 
a couple of methods, and both of the inner and outer radii are 
observable at X-ray and optical/UV bands. We chose a of 0.1 in 
order to compare the reff(r) profile and the resultant spectra with 
those in Mineshige et al. (2000). We also describe briefly the 
a-dependency using the a = 0.01, lO""* results, as well as the 
MBH-dependency comparing with 10'^(= 32)M0 BH results. 

3.1. Physical Quantities of the Flow 

Figure 1 exhibits the opacity of electron scattering (Tes), that 
of absorption (Tabs) and the effective opacity (Teff). Here, Mbh 
and a are chosen as IO^^Mq and 0.1, respectively. The fig- 
ure clearly shows that Tes (solid lines) at lOOrsch is larger 
than Tabs (dotted lines) by several orders for all the rii (1-1000). 
The deviation is much more prominent in the super-Eddington 
phases than the sub-Eddington phase {m ^ 10). Therefore, the 
local blackbody approximation must be modified, especially 
in the super-Eddington cases where Tes enormously dominates 
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over Tabs (see eq. [8]). With m ^ 100, Tgff < 1 at r^4rsch- 
Thus, the diffusion approximation for the cooling rate (eq. 6) 
will not be valid around this m. We will therefore search for 
objects with even higher rii in §4.1. The surface density S is 
proportional to t^s'- S decreases with an increasing rii as far as 
rii < 100. As m increases from 100 to 1000, E conversely in- 
creases (since increases only slightly) and the density at the 
inner region (r ^ 50rsch) is also enlarged: both of Tes and Tabs 
get larger. Consequently, the effective optical depth with m = 
1000 recovers to above unity. Open circles exhibit the radii 
where the radiative pressure (af'^/3) equals to the gas pressure 
(2kspf /mil). The former dominates over the latter within the 
inner region. The m-dependency of the radius is similar to the 
standard disk (cx ni'^/^'). 

We present the effective Compton-y parameter (eq. 10) 
for various accretion rates. Each solid line in figure 2 shows 
the distribution of as a function of radius r. The horizontal 
dotted line means of unity, above which the spectral dis- 
tortion due to Comptonization is crucial. The large increase 
of along with the increase of m arises from a decrease of 
p and an increase of f from m = 10 to 100-1000 (see fig. 5). 
Since p cx Mbh ' and T cx Mbh ''''* (as will be shown later), it is 
quite difficult to achieve a large 3'*(cx T'^-^p~^;eq. 11) with sub- 
Eddington accretion rates for all relevant Mbh- In other words, 
the huge (» 1) is obtainable only for super-Eddington cases. 




Radius (Rg^H) 

Fig. 2. — Effective Compton-v parameter, v*, for various ni. A long-daslied 
line means tliat y» = 1 , above which Comptonization is very important on the 
spectral calculations. It is clear that y, with super-Eddington accretion rate is 
much larger than the cases with sub-Eddington accretion. The inset shows the 
o-dependency of y* for m = 1000: p increases with smaller a, then y, at the 
innermost region decreases. 

We assume Keplerian rotation for the transverse Doppler ef- 
fects (§2.2). Figure 3 shows the transverse velocity is al- 
most equivalent to the Keplerian velocity vk (v^ !t 0.8 vk) up 
to the extremely large m (1000) cases. Thus, the assumption is 
not very bad. The radial velocity increases with a larger m 
(lower panel). This is the reason behind the ignition of the pho- 
ton trapping in the super-Eddington accretion. In the extreme 
super-Eddington cases {m ^ 1000), in the innermost region 
(r ^ lOrsch) is about 10% of the Keplerian one (i.e., Vi- ~ avj^). 
Thus, the inner flow could be more time variable than lower 



m cases. In the extreme limit of advection dominated flow 
(Gadv ^ 27ad)' radial velocity is expected to be ~ 0.4q;vk 
(e.g., Fukue 2000). Filled dots indicate the location of the sonic 
point, where v, equals to the sound speed c,. The sonic points 
are located around 3rsch for all m, being consistent with the 
a = 0.1 cases in Abramowicz et al. (1988). 




1 10 100 1000 10000 

Radius (Rs^h) 

Fig. 3. — Azimuthal velocity (upper) and radial velocity (lower) in the unit 
of the Keplerian velocity i'k as functions of radius. The rotation is almost Ke- 
plerian for all the in here, within ~20 % deviation at most. The radial velocity 
increases drastically with increasing m: it reaches ~10 % of the Keplerian 
velocity for m = 1000 at inner region (r ^ lOrsch)- Filled dots indicate the 
location of the sonic points (v,-/cs = 1) for each parameter set, where the sonic 
speed Cs is calculated as \J pjp. 

Next, we compare the timescales of accretion [tacc defined as 
(y- '"Sch) /v'r] with that of diffusion of photons from mid-plane to 
the surface of the disk [t^iff = /^(Tes + Tabs)/c]. Figure 4 shows 
that tacc is larger than tdiff everywhere for m < 10, while it is 
not always true for m > 100. The switch of the two timescales 
happens between m of 10 and 100. Advective cooling (photon 
trapping) begins to take place from the inner region. For exam- 
ple, the outer part (r > 40rsch for in =100 and R > 400rsch for 
m = 1000) does not suffer from the effects of photon trapping. 
Those radii can be analytically derived as 0.5 m (Begelman & 
Meier 1982). It is also clear that the accretion timescale at a 
fixed radius decreases drastically with increasing M. In super- 
Eddington cases, the structure of the inner accretion flow can 
be highly time variable. The timescale for achieving the ver- 
tical hydrostatic balance, H/cs = il^, is always much shorter 
than tacc- 

The scale height of the flow H in the unit of radius r is shown 
in figure 5 (top panel). The aspect ratio H/r increases with M, 
being geometrically thick (H/R > 0.5) at ?ri > 100. We shaU 
discuss this issue later when we are seeking ideal objects to test 
the model (§4.1). Moreover, the shape of H/R shows that the 
heating by the irradiation from the innermost part of the flow 
onto the outer radii ( ^ lO^rsch) is negligible. The radial profile 
of the mean density p is also shown (solid lines in the second 
panel). As Beloborodov (1998) noted, p at r^50rsch has a 
minimum at m ^ 100. Here, the dot-dashed line indicates a 
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Radius (Rsch) Radius (Rsch) 

Fig. 4. — Accretion timescale (tacc = (r — rsch) / '^'t', solid lines) and diffusion timescale of photons (tjiff = H{tcs +T-^bs)/c; dotted) for the slim accretion disk model. 
The viscosity parameter a is 0.1 here. The switch of the two time scale happens at m between 10 and 100. The dot-dashed lines mean the vertical, hydrodynamical 
timescale [H/cs = fl^ = 255 (A/BH/lO'M0)(f-/3rsch)'''^{''/''Sch - 1) sec for the Pseudo-Newtonian potential]. 



critical density, p(sg,z), above which self-gravity of the disk in 
the vertical (z) direction must be considered: 

p(sg,z) = ri|/(4^G) = 18(MBH/10^Mo)-2(r/3rschr' 

(r/rsch-ir'g/cml (18) 

As m increases, the relative importance of self-gravity 
[p(sg,z)/p] in the outer region is enhanced. With m > 
1000, the outer part of the flow [r^2 x lO'^rsch — 
O.OO6(Mbh/1O''-^M0)pc] will form a self-gravitating disk. Ob- 
jects with super- Eddington accretion rate are therefore good 
laboratories for studying the self-gravity of accretion flows 
by optical/Near-IR observations. Quantitative details will 
be shown in a future paper (Kawaguchi, Pierens, & Hure 
2003). Another critical density for radial self gravity, p(sg,r) 
is 16p(sg,z) (Goldreich & Lynden-Bell 1965; see Hure 1998 
for a review). The inertia of the heat (Beloborodov, Abramow- 
icz & Novikov 1997; dashed lines in the second panel) is al- 
ways much smaller than the gas rest mass density p in the 
this study (disks around non-rotating BHs): prad = af^ /c^ ~ 
10-"(r/10*'K)'*g/cm^ < p (third panel). 

Finally, the effective temperature profiles Tgf{(r) before the 
relativistic correction (solid lines) are compared with those 
with the correction (dotted lines; eq. 17). Although the flow 
with > 100 extends much closer to the central BHs than 
lower m cases, those regions are not so luminous as thought 
in the previous papers (Watarai et al. 2000, 2001; Mineshige 
2000). If the system is not face-on, the discrepancy of the 
two profiles shrinks due to the Doppler boosting and the re- 
duction of transverse Doppler effect. Viewing-angle dependent 
spectra are needed for investigation in the future. The dashed 
line also shows reff(r), but it is derived without advection (i.e., 

Gvis - 2rad)- 



reff('-) = 6.2 10- 



/ Mi 



BH 



-1/4 



V lO^Mcr 



M 



1/4 



(-) 

V '"Sch / 
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3rsch 



1/4 



K. 



(19) 
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10 100 1000 10000 

Radius (E^s^h) 

Fig. 5. — Examples of m-dependency of the physical quantities. The corre- 
sponding m are labeled. First, the scale height of the flow (from the mid-plane 
to the surface) divided by the radius is shown (top). The second panel exhibits 
the radial profile of the mean density p by solid lines. The dot-dashed line 
indicates a critical density, p(sg,z), above which self-gravity of the disk in the 
vertical (z) direction must be considered. Lower two dashed lines mean the 
inertia of the heat, Pr^j, for ni =100 and 1000. The third panel depicts the 
mean temperature f. In the bottom, the effective temperature profiles T^ffir) 
before the relativistic correction (solid lines) are compared with those with the 
correction (dotted lines; eq. 17) are shown. Dashed line also shows T^f{{r), but 
it is the case with no advection (i.e., = Q^^^)- Due to the relativistic cor- 
rection, the maximum T^(f of the slim disk is less than that of the gyj, = Q^^^ 
case. Dot-long-dashed lines in the third and forth panel present the color tem- 
perature derived with eq. 16 for the m = 1000 case. The spectral hardening 
factor can be estimated from the ratio of the color temperature to the effective 
temperature. 
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Due to the relativistic correction, the maximum T^ff of the slim 
disk becomes less than that of the Qvis = Cmd case. Dot-long- 
dashed lines in the third and forth panels present the color tem- 
perature derived with eq. 16 for the m = 1000 case. 

Since we shall use the spectra with a of 0.01 (Mbh = 
IQ^'^Mq) when testing our model with the observed spectrum, 
we here briefly discuss the a-dependency of the flow. The a- 
dependencies of all the physical quantities (figures 1-5) are 
almost the same as the standard disk: i.e., res(oc S) cx a"', 
Vr cx a, face and fjiif are proportional to a"', H/r (x a'^, p cx a"', 
T cx q;~°^^, etc. The exception arises at the innermost region, 
r < lOrsch- Namely, res(cx S), Vr and the two timescales at 
r < l.Srsch are unchanged. The transition between the two dif- 
ferent a-dependencies occurs at 1.5 < r < lOrsch- The effective 
Compton-y parameters for the three values of a (with m = 1000) 
are drawn in the inset of fig 2, showing the decrease of y^, at 
f ^ 4rsch with decreasing a. Figure 6 depicts the aspect ratio 
(top), p (second), T (third), and Teff (bottom) for different a 
(0.1, 0.01, 10"^) with Mbh = lO^'-^M© and m = 1000. The be- 
havior of T of the inner region is the opposite to that of the 
outer part. Lower a makes the self-gravity of the flow more 
important at outer radii. The sonic point approaches the center 
with decreasing a, as shown in Abramowicz et al. (1988). With 
smaller a, f{r) tends to have a hump, so that the pressure gra- 
dient [dH/dr cx d{Hf*)/dr] pushes the gas towards the central 
BH more strongly. The a-dependency of Tf^sir) is shown in 
Watarai & Mineshige (2001). 




Radius (Rsch) 

Fig. 6. — Aspect ratio (top), p (second), T (third), and T^ff (bottom) for 
different a (0.1, 0.01 and 10"^) with Mbh = lO^'^M© and m = 1000. The a- 
dependency of the physical quantities are almost the same as the standard disk, 
except for the innermost region, r ^ lOrsch- The behavior of f of the inner 
region with varying a is the opposite to that of the outer part. Arrows mean 
the evolution with increasing a. 



Computations for Mbh = 32M0 are also made in order to 
1) assess the Mbh dependency, 2) apply the present model to 
ULXs/microquasars, and 3) compare our spectra with those 
in Shimura & Takahara (1995) in the sub-Eddington regime. 
Although the ;7i-dependency of the quantities (figures 1-5) in 
the slim disk is different from the standard disk, the Mbh- 
dependency is almost identical: e.g., H/r (x Mbh'\ p(at the 
/'rad-dominated region) cx Mbh"' , T (prad region) and T^ff are 
proportional to Mbh '^"*, tacc and tdiff (piad region) are propor- 
tional to Mbh, etc. The radius (in the unit of rsch) at which 
Prad = Pgas IS proportional to Mbh^''^' ■ 

Different values of the ratio (Kabs.o/Kff.o) are used in different 
Mbh cases (§2.1), thus absorption-dominated region is located 
farther out than for Mbh = 1O^^M0. Due to the change of the 
absorption coefficient, it gets easier for the flow to have Teff less 
than unity: Teff < 1 at r < lOrsch for m = 100, and at r < 3rsch for 
m = 1000. Comparing with 10^ ^'Mq black holes, is enlarged 
by a factor of 100 at all radii and all the m. The changes of 
p and f due to the decrease of Mbh by 5 orders result only in 
0.6 order increase of (eq. 11). The rest (~ 1 .5 orders) arises 
from the decrease of (Kabs.o/'tff.o)- Since p(sg,z) cx Mbh~^, self- 
gravity of the disk around 32 M© black holes is not important 
at all. Instead, it gets more significant with larger Mbh- Again, 
Prad is always negligible. 



3.2. Emergent Spectra 

Figure 7 shows the each effect of spectral calculation on the 
emergent spectrum for the case of m = 1000 and a = 0.1. First, 
the dotted line represents the case with no advection: all the 
dissipated energy goes into radiation as in the standard accre- 
tion disks. For this curve, we put the inner radius at 3rsch- The 
long-dashed line exhibits the spectrum with advection based on 
transonic accretion flow (i.e., the same as Mineshige et al. 2000, 
except for face-on view here). We note that the local spectrum 
is still assumed as blackbody radiation here. Next, we take into 
account the relativistic correction at the innermost region of the 
flow (dot-dashed line). The flow extends inside the last stable 
orbit (3rsch; e.g., Abramowicz et al. 1988; Watarai et al. 2000; 
Mineshige et al. 2000). However, the emission from those re- 
gions, where gravity is extraordinarily strong and azimuthal 
velocity is relativistic, is drastically suppressed by the gravi- 
tational redshift and transverse Doppler redshift. As the result, 
the color temperature of the inner region (where the highest en- 
ergy photons are emitted) decreases by a factor of ^ 5.7. We 
should stress here that the bolometric luminosity of the flow L 
is not very large in spite of the super-Eddington accretion rate; 
L ~ 2.6LEdd for the relativistic corrected spectrum while it is 
5.1 Lfidd for the non-corrected spectrum. 

Next, we add the effect of electron scattering on opacity (i.e., 
the modified blackbody spectrum; short-dashed line). Finally, 
Comptonization is included in the computations (solid line). 
Radiation from the inner region is boosted towards higher en- 
ergy by a factor of ^3.4 (the ratio of the color temperature 
to the effective temperature, called as the "spectral harden- 
ing factor"; see. e.g., Shimura & Takahara 1995) in compar- 
ison with the dot-dashed line. This factor can be estimated 
if one compares the color temperature (dot-long-dashed line) 
and the effective temperature (solid line) in the bottom panel 
of fig. 5. With the effects of electron scattering (opacity and 
Comptonization), we get more gradual slopes in soft X-ray 
compared with the original spectra (dotted fine and long-dashed 
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line). The effect of Comptonization is included in the slim disks 
for the first time here (see also Shimura & Manmoto 2003 and 
Wang & Netzer 2003). 



46 



43 



1 1 1 ; 1 1 1 1 

1 oooA 


1 1 1 1 1 1 1 1 

- 








'\ \\\ 


-2 


-1 



Log(hi/[keV]) 



Fig. 7. — Each effect on the emergent spectrum for m = 1000. The dotted 
line represents the case with no advection (as in the standard accretion disks) 
with an assumed inner radius at ir^ch- Then, the long-dashed hne exhibits the 
spectrum with advection based on the transonic accretion flow (i.e., the same 
as Mineshige et al. 2000, except for face-on view here), and the bolometric lu- 
minosity of the flow L is about 5. 1 Lejij . The local spectrum is still assumed as 
blackbody radiation here. Next, we take into account the relativistic con'ection 
at the innermost region of the flow (dot-dashed line). As the result, the color 
temperature of the emission from the innermost region decreases by a factor 
of ~ 5.7, and L becomes to be 2.6 LecIcI- More, we add the effect of electron 
scattering on opacity (i.e., modified blackbody spectrum; short-dashed line). 
Radiation from the inner region is boosted towards higher energy. Finally, 
Comptonization is included in the computations (solid line). Comparing the 
solid and dot-dashed lines, the spectral hardening factor (the ratio of the color 
temperature to the effective temperature) at the innermost region is about 3.4. 

For the spectrum with all the effects (m = 1000 and a = 0.1), 
the contribution from each radial component is shown in figure 
8; r > 10^ r = 10^-100, 100-10, 10-5, 5-3 and r < 3rsch, re- 
spectively (thin solid lines). For comparison, we also draw the 
spectra without the effects of electron scattering (but with the 
relativistic correction) by thin dot-dashed lines. The spectral 
hardening factor of each spectrum is 1.3, 1.7, 2.3, 2.9, 3.4 and 
4.0 for each radial region, respectively. This means that the in- 
ner emission suffers more spectral shift by electron scattering. 
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Fig. 8.— Contributions from different radii; r > 10^ r = 10^-100, 100-10, 
10-5, 5-3 and r < 3rsch, from left to right respectively, for Mbh = lO^^M©, 
m = 1000 and a = 0. 1 . Thick soHd line shows the spectrum from the entire radii 
with all the spectral effects as same as in fig. 7, while 6 thin-solid fines mean 
the spectra from the different radial regions above. For comparison, we also 
draw the spectra without the effects of electron scattering (but with the rela- 
tivistic correction) by thin dot-dashed lines. The "spectral hardening factor" of 
each spectrum is 1.3, 1.7, 2.3, 2.9, 3.4, 4.0 for each radial region, respectively. 



Figure 9 shows the a-dependency of the emergent spectra 
with a = 0.1, 0.01 and lO'l The arrows indicate the change of 
spectra with increasing a. The meanings of two kinds of lines 
(solid and dot-dashed) are the same as fig. 7. Without the effects 
of electron scattering, the spectra are almost identical. Spec- 
tral boosting by electron scattering is quite a sensitive. Since 
Tes/Tahs fx q;^/^, an increase of a results in an increase of the 
surface temperature T (eq. 16), so we get more spectral boost- 
ing. Previous models with the effect of electron scattering in 
opacity (Szuszkiewicz et al. 1996; Wang et al. 1999) employ a 
low viscosity parameter a = . 00 1 , in order to guarantee that the 
flow is effectively optically thick. Such a low a means that the 
spectral distortion by electron scattering is very small (fig. 9). 




Log(hi/[keV]) 

Fig. 9. — Spectra with three different values of a. a = 0.1, 0.01 and 10"'. 
The arrows indicate the change of spectra with increasing a. The meanings 
of two kinds of lines (solid and dot-dashed) are the same as in fig. 7. With- 
out the effects of electron scattering, the spectra are almost identical. Spectral 
boosting by electron scattering is quite a sensitive. 

Figure 10 shows the emergent spectra with Mbh = 32Mq and 
a = 0.1. Meanings of the lines are the same as in fig. 7. Those 
spectra can correspond to those of some ULXs/microquasars. 
The spectral hardening factors are around 1,3, 1.9 and 6.5 for 
m = 1, 10 and 1000, respectively. For m = 100, it is quite dif- 
ficult to measure it, since the spectrum is highly distorted (see 
§4.2 for the reason). The m-dependency presented here is qual- 
itatively the same as the IQ^-^Mq results. Some more details 
are shown in §4.2. Due to the reduction of (Kabs.o/Kff.o) in the 
1O'^M0 calculations (§2.1), the effects of electron scattering 
upon spectral boosting are enhanced relative to 10^-^ Mq re- 
sults (fig. 7 and fig. 12). Shimura & Takahara (1995) derived 
the spectral hardening factors to be 1.7 for m = 1 and ^1.9 
for m = 10, using the radiative transfer computations includ- 
ing Comptonization. It is often assumed that the vertical dis- 
tribution of the heating rate is proportional to the volume gas 
density (Ross, Fabian & Mineshige 1992; Shimura & Takahara 
1993,1995). Our results are consistent with their values, and 
thus the method used here for spectral calculations seems to 
work quite well. Strong disk Comptonization, which is turned 
on at very high state in a couple of microquasars, is suggested 
by Kubota (2001) and Kubota, Makishima & Ebisawa (2001b). 
We note that the inferred bolometric luminosity assuming the 
standard T^^ffir) profile (e.g., Mitsuda et al. 1984) will be un- 
derestimating the true disk luminosity (perhaps by a factor of 
1.5-2), since super-Eddington accretion flows have flatter (in 
vLi, plot) spectra than the standard accretion disk (fig. 10). 
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Fig. 10. — Emergent spectra from the face-on slim disks (solid lines) with 
Mbh = 32M0 and a = 0.1. Accretion rates m used for the computations are 
1, 10, 100 and 1000 from bottom to top, respectively. The meanings of lines 
are the same as in fig. 7. Roughly, the spectral hardening factor is 1.3, 1.9 
and 6.5 for m = 1, 10 and 1000, respectively. Shimura & Takahara (1995) re- 
ported that the spectral hardening factor ~ 1.7 for m = l, and ~ 1.9 for m =10 
(Afen = 1.3 - IOM0), with a weak Mbh -dependency (^ Mbh" ')■ 

4. OBSERVATIONAL TESTS OF THE MODEL 

Here, we will test the models with the data of AGNs, for 
which detailed data at both optical/UV and X-ray bands are 
available, and some sorts of Mbh estimations have been devel- 
oped. 

4.1. The Object Ideal for The Test 

Objects ideal for the test should have the highest m, since we 
expect the largest difference in spectral shape from the spec- 
trum of the standard accretion disk. Also, the high m indi- 
cates a geometrically thick accretion flow (figure 5), then the 
self-obscuration may happen (Fukue 2000). If we observe the 
innermost region of those highest m objects (detections with 
soft X-ray will indicate it), it impUes that we are looking at 
the flow almost face-on. Thus, we can neglect the longitudi- 
nal Doppler effect (Doppler boosting), which is an annoying 
effect in the model fitting. [We simply estimate the luminos- 
ity as 47r£/L^flux, neglecting the inclination of the disk/flow. 
The inferred luminosity will overestimate the true luminosity 
somewhat in the case of a concave flow-geometry due to the 
self-irradiation (Fukue 2000; Misra & Sriram 2003).] Another 
reason to favor a large in is that the effective optical depth of the 
flow with m of 100 can be smaller than 1, so the radiative cool- 
ing rate (eq. 6) would not be accurate in that case. We should 
select even higher in objects. Another requirement for ideal test 
objects is that they should also harbor as small Mbh as possi- 
ble. Small enough Mbh allows the accretion disk to become hot 
enough for the emission from the innermost part of the flow to 
be easily observable in the soft X-ray band. 

In terms of the temporal analysis, low Mbh selection is also 
desirable in order to recognize the temporal behavior character- 
istic of the slim disk. The accretion timescale at r < 10 and < 
30 Rsch will have the order of a typical exposure time of X-ray 
satellite (^^ 10 k sec) for in= 100 and 1000, respectively (Fig- 
ure 4). Thus, the slim disk is expected to change its dynamics 
within one observation. To select the candidates for high in and 
low Mbh objects, we plot B-band luminosity of each active nu- 
clei vLy(B) as a function of Mbh as follows. 

There are several ways to estimate Mbh of AGNs. (1) the 
reverberation mapping technique provides physically mean- 
ingful Mbh estimation [MB/f(rev); e.g., Kaspi et al. 2000], 

' Mbh{[0 m]) for IC 3599 and 2E nA6+l(A(, and MsHiiev) for PG 0052+251, 



though the method is quite time consuming (an order of a year). 
Now, MBH(rev) estimations for AGNs whose FWHM(H/3) less 
than 2000 km/s have been obtained for 9 objects (Kaspi et 
al. 2000). (2) An easier way is the one assuming some 
photo-ionization models for H/3 lines [MBH(ph); Wandel, Pe- 
terson & Malkan 1999]. The two estimations agree with 
each other quite well. The method utilizes the regression line 
between the size of broad line regions /?blr and j/L^(5100 
A), determined in the reverberation mapping (Kaspi et al. 
2000): 7?BLR = 32.9[i/L^(5100A)/(10'^erg/s)]''-^ It-days. (3) 
Nelson (2000) and Ferrarese et al. (2001) found that the [O 
111] line width or the stellar velocity dispersion, cr, is corre- 
lated with MB//(rev) in AGNs, and that they overlap the bulge 
mass-to-BH mass correlation in non-active galaxies (Gebhardt 
et al. 2000; Ferrarese & Memtt 2000). Thus, the [O III] 
width also provides an estimate of Mbh, Mbh{[0 III]). We 
use the regression reported as the equation 2 in Wang & Lu 
(2001), that was derived from MsH(rev) or MsH(ph) v.s. [O III] 
width of the narrow-line and broad-line AGNs: Mbh([0 III]) 
= 10^-^'^[FW//M([OIII])/2.35/(200km/s)]3-32^Q. We con- 
firmed that an alternative regression (eq. 3) in Wang & Lu 
(2001) does not change the result significantly. 

There is some scatter around the regression lines, both in 
the Rblr v.s. h'L^iSlQOA) diagram and in the Mbh v-s. [O 
III] width plot. To minimize uncertainty, we estimate Mbh 
from the geometric mean of the three estimations, namely the 
square root of the product of MB//(rev) and MbhUO III]), or 
that of MBH(ph) and Mbh{[0 III]). In general, those evalua- 
tions agree with each other quite well: ((logjoMBH(rev or ph)- 
logioMBH([OIII]))2) = 0.54. The SDSS Early Data Release 
spectra also show a similar scatter: a ~ 0.67 in log[()MBH(ph) 
(Boroson 2003). In some objects without two available estima- 
tions, we use either of the three ' . 

For some NLSls, the [O III] line can be deconvolved into 
two components (Veron-Cetty, Veron, & Goncalves 2001, and 
references therein), namely the broad, blueshifted one and the 
narrower one with the same redshift as other lines. The latter 
seems to arise from clouds showing normal galactic rotation. 
For Mbh([0 III]), the latter [O III] width is used as Wang & Lu 
(2001) did, if the [O III] line is deconvolved into the two com- 
ponents. However, for I Zw I, it is known that neither of the 
two components appear at the galaxy redshift (Phillips 1976; 
Veron-Cetty et al. 2001): both components seem to show out- 
flows. Therefore, this object is omitted in our plot. 

In general, NLS 1 s have a narrower [O III] line width ( ^ 300 
km/s) than BLSls. To take account of a possible overestima- 
tion of the [O III] width due to the moderate spectral resolution 
(- 200km/s), Wang & Lu (2001) divide the [O III] width Hsted 
in Veron-Cetty et al. (2001) by a factor of 1.3. We follow the 
same procedure. Spectroscopic observations with higher spec- 
tral resolution and, as for I Zw I, a careful analysis on the two 
[O III] components will be required in the future to obtain a 
regression with smaller scatter. 

In addition to the objects and data in Wang & Lu (2001), 
Nelson (2000) and Kaspi et al. (2000), the following, well- 
known NLSls that are not tabulated in the literatures are 
added; namely, several NLSls with the ASCA observations (see 
Leighly 1999b). Some NLSls show giant amplitude (a factor 
of ^ 10-100) X-ray variability. We also add such objects into 
our sample (see §6.1 for more details). Table 1 lists the added 

1226+023, PG 1426+015, PG 1617+175 and PG 1700+518. 
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Fig. 1 1 . — B-band luminosity of active nuclei as a function of Mgn- Data ai'e derived from Wang & Lu (2001), Nelson (2000) and Kaspi et al. (2000). NLSls 
and narrow-line QSOs are denoted by filled circles, while BLSls and QSOs are plotted by open circles. Narrow-line objects show higher m than broad-line objects 
on average. Among the small Mbh. higher m NLSls, PG 1448-1-273 is the candidate with highest m. Then, this object can be desirable to test the overall slim disk 
model. Dashed lines are loci computed by the model for different m with a = 0.1. Optical flux has a very weak a-dependency (fig. 9). Absolute B-band magnitude 
Mb of -23 corresponds to vLv(Ji) of 10**^ erg/s. 



objects; 8 objects from Leighly (1999b) and 4 from Grupe, 
Thomas & Beuermann (2001). In total, we have 100 objects; 
74 narrow and 26 broad objects. We note that the current sam- 
ple is far from a complete sample by any means. 

The Galactic extinction is corrected using a standard extinc- 
tion curve (Savage & Mathis 1979). For all the conversion 
between optical bands and for the /^-correction, it is assumed 
that fv oc v~^-^ (e.g., Kaspi et al. 2000). Accuracy of the op- 
tical photometric data must be considered carefully: e.g., the 
nuclear fraction of the measured luminosity is 0.64-0.97 for 



i'Li,{B)> lO'^'^'^erg/s sources (Surace, Sanders & Evans 2001). 
Less luminous sources would contain less nuclear fraction. The 
optical variability is of smaller amplitude than X-ray variability. 

Figure 11 exhibits the Mbh v.s. i^Ly(B-band) plot. This 
plot will be useful to estimate the order of m (roughly cx 
Li,{B)/Mbu) in a sample of objects. The closer to the 
line of LEdd (dashed line), the higher m such objects will 
have. NLS Is/narrow-line QSOs seem to have higher m than 
BLSls/QSOs. If we do not apply the [O III]-width correc- 
tion by a factor of 1/1.3 for the Wang & Lu (2001) sample. 



Table 1 

Additional NLSls and a BLSl 



Objects Name 


z 


H/3° 


my 


4'' 

Ay 


[O III]" 




Rblk 


Mbh" 


Mbh" 






FWHM 

(km/s) 


(reL) 




FWHM 
(km/s) 




(It-day) 


(ph) 


([O III]) 


PHL 1092 


0.396 


1790 


17.0 (VVOl) 


0.13 




45.0 


168.2 


7.90 




RX J0439.7-4540 


0.224 


1010 


16.6 (G99) 


0.03 


1020 


44.7 


92.2 


7.14 


8.90 


NAB 0205-1-024 . . 


0.155 


1050 


15.4(K91) 


0.10 




44.9 


125.0 


7.30 




(Mrk 586) 
PKS 0558-504 .. 


0.137 


1250 


15.0 (R86) 


0.15 




44.9 


140.6 


7.51 




IH 0707-495 .... 


0.041 


1050 


15.7 (R86) 


0.32 


1516'' 


43.7 


18.5 


6.47 




IRAS 13224-3809 


0.067 


650' 


15.2 (Y99) 


0.23 


810 


44.3 


47.8 


6.47 


8.56 


IRAS 13349-1-2438 


0.108 


2200 


14.3 (M95) 


0.04 




45.0 


147.4 


8.02 




IRAS 20181-2244 


0.185 


370 


16.8 (ES94) 


0.24 


537 


44.5 


71.0 


6.15 


7.97 


RX J0134.2-4258 


0.238 


900 


16.0 (G99) 


0.06 


670 


45.0 


150.6 


7.25 


8.29 


WPVS 007 


0.029 


1620 


14.8 (G99) 


0.04 


320 


43.6 


17.0 


6.81 


7.23 


RX J2217.9-5941 . 


0.160 


1850 


16.2 (G99) 


0.08 


1075 


44.6 


77.0 


7.59 


8.97 


RX J2349.3-3125. 


0.135 


5210 


16.6 (G99) 


0.05 


475 


44.2 


46.0 


8.26 


7.80 



''- Line widths of H/3 and [O III] are taken from Leighly (1999b) for the objects above the horizontal line, and from Grupe et al. (1999), 
with ~ 5A~ 300 km/s resolution, for the lower objects. 

Galactic extinction Ay is derived from NED. 

log,o(i.L,(B)[erg/s]) 

log,o(MBH/M0) 

This [O III] width is uncertain due to low [O III] flux and contamination from Fe II blends (Leighly 2001, private communication). 
Then, we estimate Mbh from MBH(ph) alone. 

^ This H/3 width may be underestimated, since there seems some contamination of narrow H/3 line arises from starburst of the host 
galaxy (Leighly 2001, private communication). 

REFERENCES. — G99: Grupe et al. 1999; K91: Korista 1991; R86: Remillard et al. 1986; VVOl: Veron-Cetty & Veron 2001; Y99: 
Young et al.l999; M95: Mason et al.l995; ES94: Elizalde & Steinerl994 
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most NLSls (54 out of 74) shift to the right (larger Mbh) 
by 0.2 dex. Still there is the systematic difference of m be- 
tween NLS Is/narrow-line QSOs and BLSls/QSOs. The ob- 
jects' names are indicated for several outliers. PG 1700 and 
PG 1704 are also pointed out as the highest in objects by Collin 
et al. (2002) who estimated m of 34 AGNs from the optical lu- 
minosity andMBH(i"ev). 

The best, high in, low Mbh target is PG 1448H-273 (z = 
0.065). The inferred BH masses are MBH(ph) = 8.4x lO*" Mq 
and Mbh([0 111]) = 6.3x10^ M©. We thus estimate Mbh as 
lO*""* Mq. It has vLAB) = 10*'^* erg/s. The FWHM width of 
Uf3 and [O 111] are 820-1050 km/s and 155 km/s, respectively 
(Stirpe 1990; Veron-Cetty et al. 2001). The optical spectrum 
shows no evidence for contamination from the host galaxy into 
the optical luminosity, and it is not expected in B-band (Stirpe 
1990). These widths are still narrower than typical NLSls, in- 
dicating that this object will be one of the extreme NLSls. We 
expect the object will show the hottest accretion disk among 
the sample, which will enable us to investigate the nature of ac- 
cretion flow unprecedentedly well with soft X-ray observations 
(XMM-Newton). No X-ray observations have been performed 
for this object except for ROSAT AU Sky Survey (RASS); the 
0.1-2 keV photon index T is 3.2±0.3 (Icr), and vF^{2 keV) 
is 1.4±0.5 10-'^ erg/s/cm^ (Walter & Fink 1993). There are 
no special remarks on this object in Grupe et al. (2001) where 
some objects showing strange variabiUty during RASS obser- 
vations are reported. 

4.2. Comparison of the Models with the Available Data 

Now, we compare our model spectra (§3.2) with the currently 
available, observed data of PG 1448H-273 at optical-Soft X-ray 
bands (Figure 12). The standard disk model {in ;S 10) does not 
fit the broadband spectrum of PG 1448H-273 (dot and a tie-bow 
error), and thus a slim disk model is required. The Soft X-ray 
photon index of PG 1448H-273 (F ^ 3.2) is larger than normal 
BLSls (see e.g.. Boiler et al. 1996), but is not the largest one 
among the NLSls (F- 2.5-4.5). 

Figure 12 exhibits a comparison of spectra between the opti- 
cal (B-band) - soft X-ray (by RASS) observations and the mod- 
els. Here, Mbh is taken to be 10''-^Mq, and a = 0.1 (top) or 
a = 0.01 (bottom). The corresponding m value for each line is 
shown in the figure. The spectral hardening factors for a = 0.1 
cases are ~ 1.2, 1.3 and 3.4 for in = 1, 10 and 1000, respec- 
tively, while a = 0.01 results indicate the factor to be ~ 2.8 
and 2.3 for m = 100 and 1000, respectively. The m = 1000 and 
a = 0.1 model with all the effects (solid line) has a problem: 
the soft X-ray spectrum is too boosted towards higher energy 
to explain the observed value. We found that the PG 1448+273 
data is well fitted by a model with an extremely large accretion 
rate; m = 1000 and a = 0.01. As far as the soft X-ray spec- 
trum is concerned, both of the m= 100 and 1000 models with 
a = 0.01 can fit the data. Thus, it is very important to look at 
both of optical/UV and soft X-ray bands. We note that most of 
the dissipated energy is advected and swallowed into the BH as 
trapped photons. This is the first example where optical-soft 
X-ray spectral distribution of NLSls are explained by models. 
The dashed lines are the curves based on the previous compu- 
tation of the slim disk (Mineshige et al. 2000) with m = 1000, 
which turned out to have difficulty in explaining the soft X-ray 
luminosity/spectral slope. In the bottom panel, spectra with an 
additional hard X-ray power-law are also shown, for presenta- 
tion purpose. The photon index of 2.15 (Brandt et al. 1997; 



Leighly 1999b), with the 3000 A-2 keV spectral index of 1.6 
(Brandt, Laor & WilUs 2000) are used. 

Figure 12 also exhibits the M dependence of the model spec- 
tra with a fixed Mbh. We recall that Teff < 1 at the inner re- 
gion, with in = 100 and a = 0.1. Because of the exponential 
term in eq. 16 (the effect of the finite Teff), the surface tem- 
perature T drastically increases with the parameter set. The 
spectral hardening factor, the degree of spectral boosting, rela- 
tive to the model without electron scattering is indeed enhanced 
and becomes more strongly radius-dependent (see Shimura & 
Takahara 1995, §=3.4 for r-dependency of the factor in sub- 
Eddington AGN disks) in the m =100 case than the rii = 1000 
result. Up to m < 10, the peak frequency of the disk emission 
(Tin) varies with M in a similar fashion to that of the standard 
accretion disk: ?!„ oc M^if^l'^M^I'^ oc M^n-^l'^m^l^ Above 
that accretion rate, there is a jump in due to the ;7i-sensitive 
spectral boosting via electron scattering. Such a jump does not 
happen without electron scattering effects (dot-dashed lines). 




-2 -1 
Log(hi/[keV]) 

Fig. 12. — Spectral models of the face-on slim disk (solid lines), with 
Mbh = W^ ^Mq, for the observed data of PG 1448+273 (a dot and a tie-bow; 
errors indicate 1 ct of F and flux). Accretion rate used for the computations 
are labeled. Bolometric luminosities L are 0.021, 0.21, 1.2 and 2.6 of L^dd for 
m = 1, 10, 100 and 1000, respectively. The meanings of lines are the same 
as in fig. 7. The dot-dashed and long-dashed lines are computed without the 
effects of electron scattering: the relativistic correction is taken into account 
in the former, while not in the latter Top: The cases with a = 0. 1. Roughly, 
the spectral hardening factor is 1.2, 1.3 and 3.4 for m = 1, 10 and 1000. (For 
m = 100 it is difficult to measure it, since the spectrum is highly distorted. See 
the text for that reason.) To reproduce the optical flux, m > 1000 is required. 
The til = 1000 model with all the effects (solid line) has a problem: the soft 
X-ray spectrum is too boosted towards higher energy to explain the observed 
value. Bottom: The same as above, except for a = 0.01 here. Models for m = 
100, 1000 are shown, with the resultant spectral hardening factor of ~ 2.8 and 
2.3, respectively. For presentation purpose, we add a hard X-ray power-law 
with a photon index of 2.15 for each spectrum, with the 3000A-2 keV spec- 
tral index of 1.6. It is found that the ih = 1000 case reproduce the observed 
broad-band spectrum quite well. 

With m > 100, Tin saturates, as found in previous works 
(Wang et al. 1999; Watarai et al. 2000; Mineshige et al. 2000). 
As shown above, we have established that the slim disk 
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model is promising for PG 1448+273 for B-band and RASS 
data (m = 1000 producing 2.6xLEdd)- This is the first success 
in reproducing optical-soft X-ray emission simultaneously. To 
verify and constrain our model further, we need a much higher 
signal-to-noise ratio and a broader band pass obtainable with 
XMM-Newton. We will show the detailed comparison between 
the model and the data, which will be obtained by our ap- 
proved observation (PI; T. Kawaguchi), in a future paper. Also, 
whether we are able to fit the spectra at different flux states 
with the same Mbh provides a strong constraints on the accre- 
tion models. 

4.3. Predicted X-ray Features 

Below, we state several features predicted by the model. 
These properties will be judged by future observations. 

Is the spectrum distorted due to electron scattering? 
The most distinguishable feature of the sUm disk is the soft X- 
ray spectrum (soft excess or soft X-ray hump) from the inner- 
most region of the disk. The soft excess may not be fitted well 
by the so called multi-color disk spectra, since the shm disk 
model predicts that the soft X-ray emission is significantly dis- 
torted due to the strong radius-dependent, spectral hardening 
factor at a certain range of rh (m = 100 with a = 0.1; fig. 12 
top panel). Such a power-law-like soft excess was recently re- 
ported in a couple of NLSls, Ton S 180 (Vaughan et al. 2002) 
and Mrk 478 (MarshaU et al. 2003) with T ~ 3 and ~ 1.5 for 
soft and hard X-ray, respectively. 

Intrinsic spectral features 
The temperature of the innermost region of NLSls with 
rh ^ 1000 seems to be around 10* K. Those objects tend to 
have weaker corona than coronae of low rh objects. Then, it 
can be possible to observe the intrinsic spectral features (bound- 
free and bound-bound) of the accretion flow in the soft X-ray 
band. Such features have conomonly been attributed to "warm 
absorbers". However, the strong Comptonization (fig. 2) may 
smear those features (see Shimura 2000). Some discussions 
upon possible disk emission features in XMM data have ap- 
peared recently (Branduardi-Raymont et al. 2001; Mason et al. 
2003). 

Does the color temperature change with flux? 
The color temperature of the radiation from the innermost disk 
(Tin) may vary with flux significantly (Figure 12). Such a be- 
havior of the slim disk (significant change of Tin associated with 
less change in L) is expected to occur at rh= 10-100, due to the 
ffi-sensitive electron scattering. If we detect the significant Tin 
change, distinct from the simplest (without the effects of elec- 
tron scattering) standard disk model (Tin oc i'^"*), it is an evi- 
dence for strong electron scattering, and hence (see fig. 2) it is 
an indication of super-Eddington accretion in NLSls. 

It should be stressed that a similar Tin-change has been ob- 
served in ULXs/microquasars (see Makishima et al. 2000): e.g., 
Mizuno et al. (2001) found that fliree ULXs exhibit T]„ oc L^^, 
which is different from the relation in the standard disk. Tin oc 
Li/4. 

Is Photon trapping important or does convection dominate? 
Another expectation from the shm disk model is that the soft 
excess might be reduced at higher rh due to advection (photon 
trapping). The accretion timescale will become shorter than the 
timescale of photon diffusion (fig. 4). We may therefore ex- 
pect less soft X-ray excess-to-optical flux ratio for extremely 
high rh than for a moderately large rh case. That could be veri- 
fied from the combination of OM and EPIC observations, since 



optical/UV radiation suffers less from the photon trapping. Ab- 
sence of this effect will mean that the convective energy trans- 
port in the disk is more efficient than the radiative one. 

5. GROWTH TIMESCALE OF BHS 

The elapsed time since the gas-fueling started will have the 
order of Mbh/M (= 0.5 m~' Gy) unless a sudden change of ac- 
cretion rate occurs. Thus, NLSls and narrow-line QSOs could 
be the key objects for figuring out what turned on the efficient 
gas-accretion towards the central BHs. As described in the 
previous section, we found that some of them, including PG 
1448-1-273, have an extremely high accretion rate: rh w 1000. 
This implies that the object is really young: its age inferred 
from Mbh/M is about 10* years. Due to the rapid growth of 
Mbh, such a high rh stage will be a quite short phase. 

It may be important to note that some ULXs exhibit extended 
optical nebulae (PakuU & Mirioni 2002). Assuming that a neb- 
ula is a remnant of a supernova-like event, they estimate the age 
of one ULX as ~ 1 Myr based on the size of the nebula (~ 400 
pc in diameter), expansion velocity and H/3 lununosity. 

6. DISCUSSION 
6.1. Highly Variable AGNs 

Some NLSls show giant amphtude (a factor of ^ 10-100) 
X-ray variabihty; short-term ( ^ a day) variabihties in IRAS 
13224-3809 (BoUer et al. 1997; see also Otani et al. 1996), 
PHL 1092 (Forster & Halpern 1996; Brandt et al. 1999), 
RX J2217.9-5941 (Grupe et al. 2001), and Mrk 766 (Leighly 
et al.l996; Grupe et al. 2001), as well as long-term (~ years) 
variabilities in NGC 4051 (Guainazzi et al. 1998; Uttiey et al. 
1999), RE J1237+264 (= IC 3599; Brandt et al. 1995; Grupe 
et al. 1995a; Komossa & Bade 1999), IH 0707-495 (Leighly 
et al. 2002) and WPVS 007 (Grupe et al. 1995b), etc. To in- 
vestigate the physics and nature behind these phenomena, we 
have added such objects into figure 11; 3 NLSls, WPVS 007 
(= 1RXSJ003915.6-51 1701; Grupe etal. 1995b), RX JO 134.2- 
4258 (Grupe et al. 2000), RX J2217.9-5941 and one BLSl 
RX J2349.3-3125 (Grupe et al. 2001). RX J1304.2+0205, 
HS1702+32 (Grupe et al. 2001) are not included since tiiere 
are no optical data available. 

When rh goes up, a shorter accretion timescale is associated 
with a higher radial velocity (Figure 3 and 4). This suggests that 
such specific NLSls (and one BLSl) might have the highest rh 
among all NLSls. Some of the highly variable and transient 
AGNs (e.g., NGC 4051) may be relevant to obscuration (e.g., 
Brandt & Gallagher 2000) 

The resultant disttibution of those AGNs are shown in Fig- 
ure 13 together with other objects. It turns out that they do 
not belong to the highest th group. This is true even if we use 
one specific Mbh estimation among the three, rather than the 
geometric mean. Instead, it seems that such highly-variable 
objects are located at the intermediate region where the two 
classes merge (3 ^ m ^ 300). If this is true, the disttibution 
may indicate that such high- amplitude variability is linked with 
a ttansience between a standard disk and a slim disk due to the 
thermal instabihty (Honma et al. 1991a,b; Leighly et al. 2002). 
Further studies with a sample including more transient NLSls 
are required to examine this hypothesis. 

Janiuk, Czerny & Siemiginowska (2002) discuss the thermal 
instabihty in terms of transitions of microquasars. Transitions 
seen in some ULXs (Kubota et al. 2001a) can also be relevant 
to the instability. 
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Fig. 13. — Resultant distribution of highly variable/transient AGNs (hexagonal dots). They do not belong to the highest m group. Instead, it seems that such 
objects are located at the intermediate region where the two classes merge. A large filled square denotes NGC 1068 (§6.3). 



6.2. Photon Trapping 

As shown in fig. 4, the diffusion timescale of photons be- 
comes longer than the accretion timescale (photon trapping). 
We expect that the diffusion approximation (eq. 6) should pro- 
vide a reasonable estimation for the cooling rate, although 
Ohsuga et al. (2002) argue that it might not always be true 
for super-Eddington accretion rates. A similar trapping phe- 
nomenon can be seen for the neutrino trapping in a newly born 
neutron star (proto-neutron star) associated with a Type II (core- 
collapse) supernova. The diffusion timescale of neutrinos from 
the dense core (seconds) is longer than the dynamical timescale 
of collapse or bounce shock around the neutrino-sphere (1-10 
ms; see Burrows 1990 for a review). Numerous studies (e.g., 
Janka & Hillebrandt 1989; Mezzacappa et al. 1998, etc), adopt- 
ing the (flux-limited) diffusion approximation and Monte-Carlo 
simulations to treat the neutrino leakage, have been made for 
calculating the inner structure and light curve of the neutrino 
flux from a (convective) proto-neutron star. The long history 
of this field would help to consider the cooling rate in super- 
Eddington accretion flows. 

The inner structure of a magnetized accretion flow will be 
highly inhomogeneous (e.g., Machida et al. 2000; Kawaguchi 
et al. 2000). Photons inside the flow tend to escape by passing 
through a low-density region (see Begelman 2001). Therefore, 
the degree of photon trapping, and thus that of advection, might 
be reduced when we take this into account. However, the pho- 
ton diffusion timescale (dotted lines in fig. 4) is much longer 
than the dynamical timescale (fi^ ' dot-dashed line). Since the 
inner inhomogeneity will change within the timescale of 17^', 
photons will still have difficulty in escaping from the flow. A 
lower a enhances the deviation of the two timescales. 

6.3. NGC 1068; a Type 2 NLSl ? 

Here, we briefly discuss the accretion rate of a well-known 
Seyfert 2 galaxy, NGC 1068. The relatively narrow H/3 line 
and significant Fell multiplets in the polarized emission from 
NGC 1068 (Miller, Goodrich, Mathews 1991) lead people to 
suspect it as a type 2 counterpart of NLSl s. 

For NGC 1068, the mean [OIII] width (~1 150km/s) can not 



be used for Mbh determination, since NLR clouds seem to suf- 
fer from the interaction with the jet (Dietrich & Wagner 1998). 
The lowest FWHM (-- 200km/s) [OIII] clouds seem to be con- 
sistent with rotation around the nucleus, and the FWHM width 
of the line complex from those clouds is about 210km/s (Diet- 
rich & Wagner 1998). Consequently, the implied Mbh([0 III]) 
is 4.2 lO^'M©. Also, its polarized WI3 fine width is 3000km/s 
(Miller et al. 1991), and its intrinsic optical luminosity [i/L^(B)] 
was estimated as - 10'*'* '(/refl/0.01)-' erg/s, assuming the frac- 
tion of nuclear flux reflected into our line of sight /reo to be 0.01 
(Pier et al. 1994). Therefore, MBH(ph) = 5. 1 x 10^(/refl/0.01)-°^ 
Mq. The geometric mean is 10^-2(/refl/0.01)-°-^^ Mq. We 
note that /refl is thought to be between 0.05 (Bland-Hawthorn, 
Sokolowski, Cecil 1991) and 0.001 (Bland-Hawthorn & Voit 
1993). The inferred Mbh is consistent with the Mbh estimated 
from water masers in the central pc, (l-2)x 1O^M0 (Greenhill 
& Gwinn 1997). 

The filled square in Figure 13 indicates the location of 
NGC 1068, showing that it indeed has a high m as NLSls. The 
error-bar shown here arises from the uncertainty of /leg alone. 

6.4. Comparison with The QSO Composite Spectrum 

NLSls are the dominant source of diffuse EUV photons in 
the local universe (Edelson et al. 1999). To examine the con- 
tribution of NLSls in the distant, young universe, we need to 
understand the shape of the EUV spectrum as functions of Mbh 
andM. 

In order to recognize the difference in spectral energy distri- 
butions between NLSls and QSOs, we plot the QSO composite 
spectrum by dotted lines (Zheng et al. 1997; Laor et al. 1997; 
Telfer et al. 2002) upon our model spectra (Mbh = 1O''-^M0, 
m = 100 and 1000 with a = 0.01) in figure 14. Following Laor 
et al. (1997), the vertical normalization of the composite spec- 
trum is chosen to match with a representative NIR luminosity 
of nearby PG quasars, but re-adopting a Hubble constant of 75 
km/s/Mpc (namely we divide by 2.25 from the spectrum in fig- 
ure 6 of Laor et al. 1997). Although the m ^ 1000 period will 
be a quite short phase and thus rare, the phase with m '--^ 100 
will be common in NLSls and narrow-line QSOs (fig. 1 1). As 
shown in figure 14, the two types of objects (NLSls and QSOs) 
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show quite different distributions. The ratio of the EUV lumi- 
nosity to optical ones or EUV to X-ray ratio are really different 
in the high rii systems and QSOs. The EUV luminosity is com- 
parable in the rii = 1000 case and the QSO spectrum, while the 
EUV spectral slopes are not similar The slope in EUV band 
plays an important role in indicating to what extent EUV pho- 
tons penetrate into primordial gas clouds against self-shielding 
(e.g., Tajiri & Umemura 1998). Similar comparison can be seen 
when one plots the QSO composite spectrum onto the broad- 
band spectrum of another NLSl, Mrk 478, in Marshall et al. 
(1996). 
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Fig. 14. — Comparison between the spectral models forNLSls [solid lines; 
m =100 (lower) and 1000 (upper) with Mbh = lO^^M© and a = 0.01] and the 
QSO composite spectrum (dotted lines; Zheng et al. 1997; Laor et al. 1997; 
Telfer et al. 2002). Observational data for PG 1448+273 are also shown for a 
presentation purpose. The two types of objects show quite different distribu- 
tions. Extreme-UV luminosity is comparable in a NLSl case (m = 1000) and 
the QSO case, while optical and X-ray ones are so different. 

6.5. Other Perspectives 

As prescribed in §5, NLSls can be in the younger phase of 
central BH formation/evolution. Mathur (2000) discuss the sce- 
nario in the context of gas content and metal richness of some 
NLSls. If the "young AGN" hypothesis is the vaHd, NLSls 
provide us a good sample for studying the evolution and for- 
mation of massive BHs residing in the center of galaxies. Kro- 
ngold, Dultzin-Hacyan & Marziani (2001) compared the envi- 
ronments of NLSls with those of BLSls based on the Digi- 
tized Sky Survey, finding no statistical difference between the 
two. We will need much better spatial resolution to observe dis- 
turbances on host galaxies that may feed BHs and ignite AGN 
activity, such as a ripple structure, tidal tail, nucleus bar and so 
on. Our forthcoming observation of nearby NLSls with The 
University of Hawaii 2.2-meter telescope (PI; K. Aoki) will as- 
sess the validity of this idea. 

With recent intensive studies, it has been established that 
mass of the galactic spheroid component (galactic bulge or el- 
liptical galaxy itself) is strongly correlated with Mbh in in- 
active galaxies and AGNs (Gebhardt et al. 2000; Ferrarese & 
Merritt 2000; Nelson 2000). The tight correlation indicates 
that the growth of spheroid mass (in other words, formation 
of spheroids) and growth of Mbh (formation of BHs) are trig- 
gered by the same origin, or it implies that either controls the 
other growth/formation. The origin of the link is still unclear. 
Two, opposite, extremal ideas can be that i) Jet/Outflow due 
to gas accretion onto the central BH may govern the forma- 
tion of spheroid (Silk & Rees 1998), i.e. BHs are made before 
the bulge growth; and ii) radiation drag by massive stars dur- 
ing the bulge formation can fuel the gas towards the central 
BH (Umemura 2001), resulting in bulges formed prior to BH 



growth. Deep, optical/Near-lR, off-nuclear spectroscopic ob- 
servations of spheroids (bulges) associated with the youngest 
AGNs will reveal the origin of the close correlation. The time- 
lag between the epochs of the bulge- and BH-growth could be 
recognized, depending on the spectral features of the bulges we 
shall obtain. It will promote insights into why the two mas- 
sive, but with enormously different physical sizes, systems co- 
evolve. 

7. SUMMARY 

We have examined the effects of electron scattering (opac- 
ity and Comptonization) and the relativistic correction (gravi- 
tational redshift and transverse Doppler effect) on the emergent 
spectra from super-Eddington accretion flows (the so called 
slim disks) onto non-rotating BHs for 10'^ and 1O''-^M0 Mbh. 

The effective optical depth of the flow can be less than 
unity for rii{= M /{L^iA/c^)] = 100-1000 with 10''^Mq and 
for m = 100 with 1O''^M0, if a > O.l. An improvement in the 
knowledge of the cooling rate will be required in those param- 
eter sets. Although m-dependency and radial profile of each 
physical quantity in the super-Eddington flow are different from 
the standard disk, the MfiH-dependency of the two accretion 
modes are almost identical. The a-dependency is also similar, 
except for the innermost region (r ^ lOrsch)- narrow-line QSOs 
will be desirable laboratories for studying the self-gravity of the 
flow via optical/Near-lR observations. 

With m > 100, the spectral hardening factor via electron scat- 
tering is ~ 2.3-6.5. The color temperature of the inner most 
radiation (Tin) is not proportional to L°^^, as would be expected 
from the simplest standard accretion disk, due to the m-sensitive 
spectral hardening factor. In the face-on view, there is no jump 
of Tin if we do not take the m-sensitive electron scattering into 
account. 

The improved spectral model has been applied to the optical- 
soft X-ray emission from NLSls, for which Mbh are now esti- 
mated by a couple of ways, to evaluate their accretion rates rh. 
We have picked up one NLSl, PG 1448H-273 with an infeiTed 
Mbh of IO^ '^Mq, among the highest m candidates. The broad- 
band spectral distribution was successfully reproduced for the 
first time by the model with m = 1000 and a = 0.01. 

This extremely high m implies that this object, as well as 
some other highest m systems, is really young: its inferred age 
(Mbh/M) is less than 10^ years. 

We have also briefly discussed the accretion rate m of tran- 
sient and highly variable NLSls, finding that those are located 
at 3 ^ m ^ 300. Such a moderately high accretion rate can in- 
dicate that this variability is relevant to the thermal instabiUty. 
A possible type 2 counterpart of NLSls, NGC 1068, seems to 
have m similar to NLSls. 

A more accurate estimate for absorption opacity, for which 
we assume 30 times free-free opacity in this study for simplic- 
ity, is required for better spectral calculations. It is because the 
absorption opacity determines the location of the last thermal- 
ization surface, and hence the amount of redistributed photons 
by Comptonization. Also, viewing-angle dependent spectral 
models are needed. 

This paper could never have been produced in this way with- 
out the use of the numerical code for solving the radial struc- 
ture of super- and sub-Eddington accretion flows/disks that has 
been constructed and developed by Ryoji Matsumoto, Fumio 
Honma and Mitsuru Takeuchi. The author wish to thank Karen 
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Toshiya Shimura, Agata Rozanska, Jerzy Madej, Bozena Cz- 
emy, Ed Baron, Jean-Marc Hure, Suzy CoUin, Chris Done and 
Omer Blaes for helpful discussions. He also thanks Brandon 
Carter for patient and careful reading of the draft, Shin Mi- 
neshige for continuous encouragement, and an anonymous ref- 
eree for careful reading and helpful comments. The main com- 
putations were performed at the Department of Physics and 
Astronomy of the University of Oklahoma, when the author 
was there as a visiting graduate student. He appreciates all the 



member of the Department, especially C. Matsumoto and K. 
Leighly, for their hospitality. This research has made use of 
the NASA/IPAC Extragalactic Database (NED) which is oper- 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Administration. This work was supported in part by the 
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